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SUNir^RY 


Laser-Doppler measurements are reported for laminar and turbulent flow 
through a ninety-degree bend of circular cross section with mean radius of 
curvature equal to 2.8 times the diameter. The measurements were made in cross- 
stream planes 0.58 diameters upstream of the bend inlet plane, in 30, 60, and 
75 degree planes in the bend and in planes one and six diameters downstream of 
the exit plane. Three sets of data were obtained: for laminar flow at Reynolds 

numbers of 500 and 1093 and for turbulent flow at the maximum obtainable Reynolds 
number of 43000. The results show the development of strong pressure driven 
secondary flows in the form of a pair of counter rotating vortices in the 
streamwise direction. The strength and character of the secondary flows were 
found to depend on the thickness and nature of the inlet boundary layers: inlet 
conditions which could not be varied independently of Reynolds number. The 
quantitative anemometer measurements are supported by flow visualisation studies. 
Refractive index matching at the fluid-wall interface was not used; the 
measurements consist therefore of streamwise components of mean and fluctuating 
velocities only supplemented by wall pressure measurements for the turbulent 
flow. The displacement of the laser measurement volume due to refraction is 
allowed for in simple geometrical calculations. The results are intended for 
use as benchmark data for calibrating flow calculation methods. 


This research was supported by contracts from the CEGB Central Electricity Research 
Laboratories, London, England, and from the NASA Lewis Research Center, Cleveland, 
Ohio, U.S.A. 
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1 . INTRODUCTION 


As a fluid flows through a straight pipe or duct and into a bend the 
pressure, which in the straight section is uniform across the flow, must adjust 
in the bend to counter centrifugal forces. The pressure is greatest at the 
outer wall farthest from the centre of curvature and least at the inner wall 
nearest the centre of curvature and the cross-stream pressure gradient in the 
bend has two well known effects on the flow. At the bend inlet the boundary 
layer on the outer wall experiences the effects of a positive streamwise 
pressure gradient which may in a tight bend be strong enough to produce local 
separation; conversely, the inner wall boundary layer is accelerated. The 
reverse occurs at the pipe exit where local pressure gradients of the opposite 
sign appear as the flow adjusts to uniform pressure conditions downstream. 
Inside the bend the presence of side wall boundary layers ensures that the 
conditions for radial equilibrium cannot be satisfied exactly over the entire 
field. The slow moving fluid in these layers is forced inwards towards the 
centre of curvature and continuity ensures that the high speed fluid near the 
axis is driven outwards. The result is a secondary flow superimposed on the 
main flow. 

If the bend is of circular cross section, or its height is not too 
much greater than its width measured in the radial direction, the effect of 
the secondary flow is to displace the region of maximum velocity from the 
centre towards the outer wall. If, however, the boundary layers are thin and 
the secondary flow is weak the velocity distribution in the inviscid central 
core will approximate to that of a free vortex with velocity increasing 
towards the centre of curvature. The effects of longitudinal curvature on 
the boundary layers on the inner and outer walls also tend to counteract the 
effects of pressure-dri ven secondary flow. In the inner boundary layer the 
shear stress is reduced by curvature and the boundary layer grows relatively 
slowly. In the outer boundary layer the shear stress is increased and the 
layer is retarded and thickened. Thus in the flow in a high narrow curved 
channel, where the secondary flows are weak and confined to the edge zones, 
the region of maximum velocity is found closer to the inner wall (1). The 
effects of curvature on laminar boundary layers are weak in ducts of moderate 
aspect ratio and are masked by the secondary flow; the fractional change in 
shear stress is of order 5/R (where 6 = boundary layer thickness and 
R s radius of curvature) . The effects of curvature on the turbulence in 
turbulent boundary layer flow are an order of magnitude greater (2) and may 
be large enough to counteract the effects of secondary flow and to play a 
significant role in determining the flow structure. 

The presence and origin of secondary flow in bends have long been 
known and understood although obvious measurement difficulties associated with 
probe explorations have hindered detailed investigation and have largely 
restricted experiments to flow visualization studies and overall pressure loss 
measurements. The former (3) show that a laminar boundary layer in a ninety- 
degree pipe bend is deflected through about 45 degrees while the deflexion 
of a turbulent layer does not exceed 30 degrees as might be inferred from the 
different inlet profiles. Work done prior to 1970, which consists almost 
entirely of wall pressure measurements and pressure probe surveys of the outlet 
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plane, is summarized and reviewed by Ward-Smith (4). These data provide a 
tn h! ^^0*^ /'JPT'ical pressure drop correlations but are necessarily too limited 
to be used for the validation of flow calculation methods. The development 

standard design tools 

has, fortunately, been paralleled by the development of the laser Doppler 

ol^viating the need for probe insertion, has made possible 
wre detailed measurement of complex flows including those in bends. Humphrey 
lami^ measurements of the longitudinal velocity component in 

is extendJd^-i^rfil^t a strongly curved duct of square cross section. This work 

™a„ ill 

® P^’P®» Which is the subject of the 

present contribution, presents additional measurement difficulties and has 
consequently received less attention. Rowe (8) reports total -pres sure 

aoolal^^trbP laser-Doppler measurements 

appear to be limited to the laminar flow experiment described by Agrawal et 

are rifraJtPH^irthp because in off-diameter planes the laser beams 

refrlcnv^i^Hi^Pc^ !f"less care is taken to match the 

indices. If this is not done the measurements will be limited to 

ioW i^Hn^ni^rpII^II^J fluctuating velocity. The measurement 

volume IS displaced but to a position which can be located quite accurately 
by simple calculation (10). au^-uraueiy 


finui P^P®^ we report measurements of laminar and turbulent 

flow in and downstream of a tight 90-degree pipe bend. The results which 

f < = • ^ usefSl coninEiMon 

to the data bank available for the further testing and development of 
numerical flow calculation methods. 


sions . 


Professor J.H. Whitelaw contributed valuable advice in numerous discus- 


2.1 The test bend 


2 . EXPERIMENTAL ARRANGEMENT 


The geometry and dimensions of the test bend are shown in Figure 1 The 
dimensions are similar to those of the square-sectioned bend of (6.7) with 
radius of curvature equal to 2.8 times the internal diameter. The bend was 
machined accurately from two halves of a rectangular plexiglass block split 
^ e plane of synmetry. It was thus possible to maintain a tolerance of 

extern^ diameter and to minimise refraction problems at the 

, air plexiglass, interface by the retention of a rectangular external 
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section. The bend was fixed in the horizontal plane downst^am of a straight 
tube of the same diameter (48 mm) and 240 mm long; a second straight tube, 
480 mm long, was fitted downstream. Photographs of the overall arrangement 

are shown in Figure 2. 

Figure 3 shows diagrammatical ly the closed circuit water 
which the bend was fitted immediately downstream of a contraction nUed to 
the plenum chamber. The flow rate was constantly monitored using Fischer and 
Porter precision flowmeters. The maximum obtainable Reynolds number, — — ^ » 
was 43000. For the turbulent flow measurements the water was seeded 
with minute quantities of milk in order to increase the scattering particle 
concentration and, consequently, the particle arrival rate (11). 


2.2 Instrumentation and signal processing 

The optical arrangement, which is similar to that described by Taylor 
et al (7), is shown diagrammatically in Figure 4. The light beam generated 
by a 5 mW helium-neon laser is focussed on to a radial diffraction grating, 
collimated in order to reduce spherical aberration and then focused to fom 
the scattering volume at the point of measurement in the bend. Measurements 
in the vertical plane were made in the straight pipe sections at each end of 
the bend. For these the light beam was turned through 90 deg using a mirror 
which was optically flat to within one tenth of the wavelength of light. 

Forward scattered light collected by the objective lens is focussed on 
passes through a pinhole to an EMI 9658 photomultiplier, the output from which 
is processed in a Cambridge Consultants CCOl frequency tracking demodulator 
and associated equipment shown in the block diagram of Figure 5 The principal 
characteristics of the laser-Doppler anemometer are given in Table 1 below. 


Table 1 

Ootical Characteristics of the Laser Doppler 

Anemometer 

» 

Focal length of imaging lens 

300 mm 

Half-angle of beam intersection 

9.11° 

Fringe separation (line-pair spacing) 

1 .9 pm 

Number of fringes in measuring volume 
Intersection volume diameter calculated at 

84 

1/e^ intensity 

Intersection volume length calculated at 

0.167 mm 

1/e^ intensity 

1.391 mm 

Photomultiplier pinhole diameter 

0.50 mm 

Transfer constant 

0.500 MHz/ms 
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2.3 Scope and accuracy of the measurements 

Refraction of the laser beam at the water-plexiglass interface not only 
displaces the measurement volume by an amount which can be readily calculated 
but also, and more importantly, prevents the measurement of cross-section 
velocity components. The method employed successfully for the latter in the 
rectangular bend experiments (6, 7) fails for the round pipe because, unless 
the refractive indices of the fluid and pipe wall are equal, the twin laser 
beams entering the pipe in a plane at 45 degrees to the axis cannot converge. 
Consideration has been given to the use of fluids other than water which would 
allow the indices to be equalized. Unfortunately it turns out that such fluids 
either possess highly unpleasant or toxic properties or are so viscous that it 
is difficult to realize Reynolds numbers high enough for turbulent flow. A 
case in point is the experiment described by Agrawal et al (9) where the use 
of a glycerine-water mixture to match the indices limited the investigation to 
laminar flow. Since the main objective of the present work is the study of 
turbulent bend flow, the limitation on the scope of the measurements to stream- 
wise velocities only has had to be accepted at this stage. 

There is no difficulty in measuring the streamwise velocities or in 
calculating the displaced position of the measurement volume. The measurements 
are made at points along the optical axis which is refracted at the pipe wall 
at an angle to the plane of symmetry to which it was originally parallel. A 
straightforward calculation produces the beam trajectories shown in Figure 6. 
The actual location of the measurement points along the optical axis is 
obtained by applying the correction described by Boadway and Karahan (10) for 
the effects of refraction due to the bend curvatures. These two corrections 
have been applied independently to the results presented here. It was also 
found that signal quality was only slightly impaired by refraction and that 
measurements were feasible for points in the flow distant more than one tenth 
of the radius from the wall in the vertical plane. 

Estimated measurement errors are summarized in Table 2. Corrections 
were not applied for the effects of finite transit time and instrument noise 
broadening. Systematic errors of up to 2.5 per cent are mainly associated with 
velocity gradient effects near the walls and have been included in the table. 


Measurement errors 

Table 2 


Quantity 

Systematic error 

Random error 

X 

±0.5 mm 

±0.02 mm 

y.z 

±0.2 mm 

±0.02 mm 

6 

±0.3° 

±0.170 

6 

±0.05° 

nil 

Ug (laminar) 

±0.8% 

±0.8% 

Ug (turbulent) 

±1.5% 

±0.5% 

U/Ug 

up to 2.5% 

±1.5% 

u'/Ub 

up to 3% 

up to ±3% 


5 



The measurements were made in cross-stream planes 0.58 diameters 
upstream of the bend inlet plane, in 30, 60, and 75 degree planes in the bend 
and in planes at one and six diameters downstream of the exit plane. Three sets 
of data were obtained: for laminar flow at Reynolds numbers of 500 and 1093 

and for turbulent flow at the maximum obtainable Reynolds number of 43000. 


3. RESULTS AND DISCUSSION 


3.1 Flow visualization 

Observation of hydrogen bubble paths in low speed flow revealed no 
unusual features such as the local separation reported in (5) for fully 
developed flow through the square section bend. The visual observations were 
entirely consistent with large scale secondary flow patterns detected by 
laser Doppler anemometry. Hydrogen bubbles were generated in the flow by a 
thin wire electrode inserted through static pressure tappings, the trajectories 
were photographed with typical results shown in Figure 7. Figures 7.1 and 7.2 
show bubble paths at the 15 and 30 degree stations where the inward movement 
of low velocity fluid in the upper part of the bend is clearly visible. The 
secondary flow appears to be directed inward at a greater angle at the 15 
degree station than further downstream at 30 degrees. The corresponding outward 
movement of the core flow is shown by the distinct trace from the electrode 
tip. Figure 7.3 shows the hydrogen bubble sheet generated near the outer wall 
at the 30 degree station; the flow here is upward and along the wall toward 
the inside of the bend. 

3.2 Laminar flow measurements 

The results of velocity measurements at a Reynolds number of 500, mean 
velocity 10.5 mm s"1 , are presented graphically in Figures 8 and 9. The data 
are tabulated in Appendix 1. Figure 8 shows the inlet conditions: horizontal 

and vertical velocity profiles 0.58 diameters upstream of the bend. The 
boundary layer thickness is about 0.27 d so that there is a significant 
potential core. The apparent slight asyirmetry shown in the results is withm 
the estimated error band. Contours of the streamwise velocity in four cross- 
stream planes are presented in Figure 9 and show the effects of strong 
secondary flow which displaces the region of maximum velocity to the outside 
of the bend. The secondary flow develops gradually: in the 30 degree plane 

the only visible effect is a thickening of the shear layer on the inside of 
the bend. By 60 degrees the secondary flow has become almost fully developed 
as is shown by the distortion of the contours. The fourth set of contours, 
in the plane one diameter downstream, shows that the secondary motion persists 
without significant decay to this distance. 
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sffects of doublInQ the mesn velocity dnd Reynolds number to 
23 mm s~' and 1093 respectively are shown in Figures 10 and 11. The data are 
tabulated in Appendix 2, The inlet boundary layers are reduced in thickness 
by about fifteen per cent to 0.23 d, the transverse pressure gradient is 
increased and the secondary flow is intensified. The peak velocity in the 30 
degree plane is displaced from the plane of symmetry but remains closer to the 
centreline. This displacement, which persists downstream, has also been 
observed in the flow in a square bend (7). The shape of the contours near the 
outside of the bend reflect scatter in the data attributable mainly to 
difficulties in locating accurately the measurement volume near the wall and 
partly, perhaps, to slight asymmetry in the flow. 

In neither case do the measurements show any unexpected feature: 
the secondary flow superposed on the main motion is essentially as described 
qualitatively in references (1) and (3); it develops gradually and persists 
downstream after the removal of the transverse pressure gradient responsible 
for its generation in the bend. The strength and location of the pressure 
driven secondary flow depend on the inlet conditions and the Reynolds number 
which could not be varied independently in the present experiments. The 
stronger secondary flow observed at Re = 1093 must be attributable mainly to 
the increased centrifugal forces and transverse pressure gradient rather than 
to the relatively small changes noted in the initial conditions. The present 
contribution is of quantitative velocity data in the bend: these are 

considered to be accurate enough to be useful in checking the numerical accuracy 
of flow prediction methods. 


3.3 Turbulent flow measurements 


At the maximum flow velocity of 0.92 m s"!, Reynolds number 43000, 
the inlet boundary layers are turbulent. The static pressure variation on the 
bend wall is large enough to be measured accurately with the results shown in 
Figure 12. At the bend inlet the negative pressure gradient on the inner wall 
IS roughly twice the positive gradient on the outside. These initial gradients 
resulting from the change from straight to curved flow disappear at about 25 
degrees so that in the mid-section of the bend (250<6 <75°) a quasi equlibrium 
condition is reached with approximately uniform pressure on the inner and 
outer walls. The transition to straight flow downstream is signalled by the 
appearance of strong gradients for 6 > 75 deg. The overall pressure drop 
IS about 0.3pU|. ^ 

Laser Doppler measurements of the mean and r.m.s. fluctuation of the 
streamwise velocity at inlet are shown in Figure 13 and contours of these 
quantities at downstream stations in Figures 14 and 15. Numerical values are 
tabulated in Appendix 3, The inlet boundary layers of Figure 13 are much 
thinner, 0.09d, than in the two laminar flows investigated and the presence of 
a large central region of uniform velocity flow influences significantly the 
development of secondary flow downstream. The mean velocity contours presented 
in Figure 14 show that by the 30 degree plane the potential core has adjusted 
to the cross-stream pressure gradient imposed by the bend and the region of 
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maximum velocity has been displaced toward the inside. The results also 
differ from the corresponding laminar flow measurements in that the thickness 
of the outer wall boundary layer is increased relative to that on the inside 
of the bend due mainly to the effects of the potential flow adjustment and the 
wall pressure gradients of opposite sign but also, quite possibly, to the 
effects of longitudinal curvature on the turbulence structure. In the 60 degree 
plane distortion of the contours suggests the presence of strong secondary 
flow largely confined to the region of steep velocity gradients near the inner 
wall. This secondary flow grows with distance through the bend and by the 75 
degree plane extends over and beyond the inner half of the section. The removal 
of the transverse pressure gradient in the straight flow downstream results in 
a shift of the velocity maximum toward the outside. A large secondary flow, 
comparable with those observed in the laminar flow experiments, now extends 
over the entire flow area and persists in the plane six diameters downstream of 
the bend where, however, the transverse velocity gradients have decayed to about 
one half of their values one diameter downstream. 

The contours of turbulence intensity presented in Figure 15 are 
consistent with those of the mean velocity. As is to be expected high values 
of the velocity fluctuation are associated with steep mean velocity gradients 
and with the distortion of the symmetrical inlet flow. 

These results may be compared with the measurements made by Humphrey, 
Whitelaw andYee(6) in a tight bend of square cross-section. The inlet 
boundary layers of that flow were substantially thicker than in the present 
flow and contained stress driven secondary flows towards the corners. As a 
result the pressure driven secondary flow in the bend extended over the entire 
flow area as in the two laminar flows of the present study. 


4. CLOSURE 

The measurements show, as expected, the presence in bend flow of 
strong cross-stream flows produced by the inward motion of slow moving fluid 
in a radial pressure gradient. The strength of the secondary flow is influenced 
by the initial inlet flow conditions: when the inlet boundary layers are thin, 

as in the turbulent flow case reported here, the initial flow adjustment is an 
inward shift of the velocity peak in the potential core. The secondary flow 
then develops initially in the inside of the bend. For all cases the 
secondary flow takes some time to develop and its effects are still increasing 
by the exit and persist downstream. 

The effects of refraction at the doubly curved fluid-wall interface 
prevent the measurement of cross-stream velocity components. The investigation 
of flows in the round bend to the same level of detail as has been achieved 
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previously in square bends will require refractive index matching. For the 
streamwise velocity measurements reported here simple geometrical corrections 
allow for the displacement of the measurement volume by refraction. 

/laminar flow data may be used to check the numerical accuracy of 
f ow prediction methods; the turbulent flow measurements provide the additional 
information needed to validato mathematical models of turbulence. 
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Tabulated Data for Reynolds Number 500 Measurements 

All velocities normalised with the bulk velocity, Ug = 0.0105 m/s. 
Symbols as defined in Figures 1 and 6. 


x/d 

y/r 

z/r 

U/Ug 

y/r 

z/r 1 

U/Ug 

-0.58 

0.9 

0 

0.509 

0 

-0.9 

0.539 


0.8 


0.865 


-0.8 : 

0.820 


0.7 


1.167 


-0.7 : 

1 .106 


0.6 


1 .393 


-0.6 ^ 

1.335 


0.5 


1.546 


-0.5 i 

1 .484 


0.4 


1 .603 


-0.4 i 

1 .595 


0.3 


1 .648 


-0.3 

1 .655 


0.2 


1.66 


-0,2 

1.673 


0.1 


1 .664 


-0.1 

1.681 


0.0 


1.663 


0.0 

1 .687 


-0.1 


1 .658 


0.1 1 

1.682 


-0.2 


1.66 


0.2 I 

1 .678 


-0.3 


1.648 


0.3 

1.667 


-0.4 


1.607 


0.4 

1 .639 


-0.5 


1.55 


0.5 

1.552 


-0.6 


1.408 


0.6 

1.419 


-0.7 


1 .18 


j 0.7 

1.167 


-0.8 


0.873 


0.8 

0.847 


-0.9 


0.505 


0.9 

' 0.482 

1 
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Appendix 2 

Tabulated Data for Reynolds Number 1093 Measurements 


All velocities normalised with the bulk velocity, Ug = 0.023 m/s. 
Symbols as defined in Figures 1 and 6. 
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APPENDIX 4 


DEFINITION OF SYMBOLS 


Roman Characters 

Cp Pressure coefficient (figure 12) 

d Pipe diameter 

lii Mass flow rate 

p Static pressure 

*^ref Static pressure one diameter upstream of pipe bend 

R Radius of curvature of the bend 

r Internal radius of pipe = -^ d 

Re Reynolds number Ugd/v 

U Streamwise mean velocity 

Ub Bulk velocity 4m/(7rpd^) 

R.m.s. streamwise velocity fluctuation 
X Distance measured along the pipe and bend centreline 

y Distance measured horizontally from centreline (figure 1) 

z Distance measured vertically from centreline (figure 1) 

Greek Characters 

3 Half angle of beam intersection 

6 Boundary layer thickness 

9 Bend angle (figure 1) 

<t> Angle substended at pipe centreline (figure 12) 

V Ki nemati c vi scos i ty 

P Density 
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Figure 3 Closed circuit water flow rig (plan view) 







Figure 5 Block diagram of Doppler signal processing system 
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Figure 6(ii) Cross-section through upstream and downstream tangents 





Figure 7.1 Flow visualization: Re = 1093. Flow from top left 

to top right; wire inserted at 6 = 15°, y/r = 0, z/r = 1.0 



Figure 7.2 Flow visualization: Re = 1093. Flow direction from top 

left to top right; wire inserted at 0 = 30°, y/r = 0, 
z/r = 1 .0 
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Figure 7.3 Flow visualization: Re =1093. Flow direction from top 
right to top left, wire inserted at 0 = 30 , y/r = 1.0, 
z/r = 0 



Figure 8 Horizontal, 0,and vertical, A, profiles of mean velocity 
in the laminar approach flow 0.58 diameters upstream 
of the inlet plane. Re = 500 
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Figure 13 Profiles of mean velocity U/Ug and turbulence intensity 

u'/Un in the turbulent approach flow 0.58 diameters upstream 
of tne bend. O horizontal , A vertical . Re = 43 000 
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